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Abstract 

In this paper the characteristic matrix method is used to study the propagation of electromagnetic waves 
through one-dimensional lossy photonic crystals composed of negative and positive refractive index material 
layers with symmetric and asymmetric geometric structures with a defect layer at the center of the structure. 
First, the positive index material defect layer is considered, and the effects of the polarization and the angle of 
incidence on the defect mode in the transmission spectra of the both structures are investigated. The results 
show that the number of the defect modes in the transmission spectra depends on the geometry (symmetric 
or asymmetric) of the structure. In addition, it is shown that the defect mode frequency increases as the 
angle of incidence increases. This property is independent of the geometry of the structure. Then, for normal 
incidence, the negative index material defect layer is considered, and the properties of defect modes for both 
structures are investigated. The results can lead to designing new types of transmission narrow filters. 
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I. INTRODUCTION 



Photonic crystals (PCs) are artificial dielectrics or metallic structures in which the refractive 
index changes periodically. Due to their unique electromagnetic properties and important scientific 
and engineering applications, research on PCs has been an attractive topic in optics during past 
two decades [lHg]. The spatial periodicity of the PCs causes a range of forbidden frequencies 
called photonic band gap [7|, |_Sf| . The gap is formed as a result of Bragg scattering in a periodical 
dielectric structure through which no electromagnetic wave can propagate. By introducing a layer 
with different optical properties which breaks the periodicity of the PCs structure, localized defect 
modes will appear inside the band gap, very similar to the defect states that are generated in the 
forbidden band of doped semiconductors. The appearance of the defect modes is due to the change 
of the interference behavior of electromagnetic waves. 

In 1967 Veslago j{J] predicted the existence of materials with negative refractive index called 
metamaterials or Negative Index Materials (NIMs). NIM also known as left-handed material, has 
negative permittivity and permeability simultaneously. After experimental realization of metama- 
terials by Smith et al. [h]], such materials have received extensive attention for their very unusual 
electromagnetic properties fill 17]. Recently, with the possibility of producing metamaterials, PCs 
with metamaterials, called metamaterial photonic crystals (MPCs) have been made. In 2003, Li 
et al. have interestingly reported the appearance of an additional gap called zero-re gap in the 
transmission spectra of a one-dimensional (ID) PC composed of NIM and positive index material 
(PIM) dielectric layers. 

In several papers the properties of the defect modes in different ID conventional PC and ID 
MPC structures are reported [19H34I]. Following the interesting report by Wu et al. [24| on the 



properties of the defect modes in ID conventional PC with symmetric and asymmetric structures, 
here we investigate the properties of the defect modes in the transmission spectra of ID symmetric 
and asymmetric lossy MPCs structures with the defect layer at the center of the structure. First, 
we consider the PIM defect layer, and after that, the NIM defect layer has been considered. 

The outline of this paper is as follows. In Section 2, two PMC structures, the permittivity and 
permeability of the NIM layer, and also the characteristic matrix method and its formulation are 
presented. The numerical computation and the results are given in Section 3, and the paper is 
concluded in Section 4. 
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II. MPC STRUCTURE AND CHARACTERISTIC MATRIX METHOD 



A ID MPC with asymmetric and symmetric structures in air with a defect layer at the center 
of the structure are shown in Figure la and Figure lb, respectively. A is a NIM layer, B is a 
PIM layer, and defect layer C is either PIM or NIM. The NIM is assumed to be dispersive and 
dissipative. The permittivity and the permeability of layer A are complex, and are given by 



saU) = 1 + 



+ 



0.9 2 - f 2 - 17/ 11.5 2 - f 2 - i-yf 



(1) 



Ma(/) = 1 + 



(2) 



0.902 2 - p - iff 

where / and 7 (given in GHz) are respectively frequency and damping frequency. Plots of the real 
parts of the permittivity and permeability of layer A, e' A and fi' A , versus frequency for four different 
loss factors are shown in Figure 2. Also, more details of the various aspects e' A and \1 A have been 



discussed in our previous study 



35]. 



36[|, which is the most effective 



Our calculation is based on the characteristic matrix method 
method to analyze the transmission properties of PCs. In the absence of the defect layer, the 
characteristic matrix of the periodic structure (AB) is given by: 
M[d] = (M A M B ) N 

where N is the number of the lattice periods. In the presence of the defect layer, the characteristic 
matrix for asymmetric (AB) N I 2 C{AB) N / 2 , and symmetric (AB) N / 2 C(BA) N I 2 structures, respec- 
tively, are: 

M[d] = (M A M B ) N / 2 M C (M A M B ) N / 2 , and M[d] = {M A M B ) N I 2 M C {M B M A ) N I 2 , 

where M A , M B and Mq are the characteristic matrices of layers A, B, and C. The characteristic 



matrix M, for TE waves at incidence angle 9q in vacuum is given by 



Mi 



cos 7 f — sin 7; 



(3) 



-1 Pi sni7j cos7i 

where 7.; = (w/c) njdjcos^j, c is speed of light in vacuum, 9i is the angle of refraction inside the 



layer i with refractive index m and pi = yjsil \i{ cos#j, where cos#j = yl — (n 2 , sin 2 #o/ ra i)> i n 
which no is the refractive index of the environment where the incidence wave tends to enter the 
structure. The refractive index is given as rij = 3l], [37 j , where the positive and the negative 

signs are assigned for the PIM and NIM layers, respectively. The final characteristic matrix for an 
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N period structure is given by: 

[M(d)f = J] M, = [ ) (4) 




i=l,...,iV 

where niij(i,j = 1,2) are the matrix elements of [M(d)] 7V . The transmission coefficient of the 
multilayer is calculated by: 

t = 2 J^ ( 5 ) 

(mn + mn Ps) Po + (m 2 i + m 22 p s ) 

where po = uq cos 9q an d p s — n s cos 9 S , with n s being the refractive index of the environment 
where the wave leaves the crystal with angle S . The transmissivity of the multilayer is given by 

T = ^|t| 2 . (6) 
Po 

The transmissivity of the multilayer for TM waves can be obtained by using these previous 
expressions with pi = cos 0j, p = £ ^ a , and p s = 



III. NUMERICAL RESULTS AND DISCUSSION 

Based on the theoretical model described in the previous section, the transmission spectrum 
of the lossy MPC structures is calculated. In the first part, we use equations (1) and (2) for the 
permittivity and the permeability of layer A. Layers B and C are assumed to be PIM, and the 
refraction indices are ng = 1 and nc = 2.7. The thickness of layers A, B, and C are considered as 
d>A = 6 mm, ds = 12 mm, and dc = 24 mm. Calculations a re p erformed with N = 16 



38, 



39]. 



351 ] in the 



s k 

The transmission spectra of TE and TM polarized waves for the asymmetric structure at various 
angles of incidence and for 7 = 0.2 x 10~ 3 GHz are shown in Figures 3 and 4, respectively. As it 
is seen, the defect layer causes two different modes. The modes are identified by numbers 1 and 2 

n 

in the figures. This is in sharp contrast to the report by Wu et al. [241 ] for the defect modes in ID 
PCs with PIMs where only one defect mode for the asymmetric structure has been observed. 

In this asymmetric ID MPC, the frequency dependence of the defect modes as a function of the 
incidence angles for TE and TM waves and for 7 = 0.2 x 10~ 3 GHz is illustrated in Figure 5. TE 
and TM modes start from the same value for zero incidence angle, as expected, but the starting 
points are different for different mode numbers. As the incidence angle increases, the frequency 
of the TM defect modes increases, but the frequency of the TE modes remains almost unchanged 



24j in the study of the 



for both defect modes. Similar behavior has been reported by Wu et al. 
symmetric PIM geometry structure. It is interesting to mention that the first TM defect mode 
disappears for an incidence angle greater than ~ 50°. 

Now, the properties of the defect mode in the transmission spectrum of the ID MPC with 
symmetric structure with a PIM defect layer at the center are investigated. The results for TE 
and TM polarized waves for different incidence angles and for 7 = 0.2 x 10~ 3 GHz are shown in 
Figures 6 and 7, respectively. As it is seen from the figures, there exists only one defect mode in 
accordance with the work reported by Wu et al. [2J] in the asymmetric structure. 



The frequency dependence of the defect mode on the incidence angle for TE and TM waves is 
shown in Figure 8. It is seen that the frequency of the defect modes increase with a higher rate for 
TM wave were the incidence angle passes over ~ 35°. 

In the last part, properties of the defect mode in the transmission spectra of both geometric 
structures with an NIM defect layer at the center and for normal incidence case are investigated. 
The asymmetric and symmetric structures which were used before are modified by replacing an 
NIM defect layer (layer C) in which the permittivity and permeability follow the equations (1) 
and (2), exactly like layer A. The other parameters are kept the same as the first part, PIM 
defect layer. The results of the asymmetric and symmetric structures for 7 = 0.2 x 10 -3 GHz are 
respectively shown in Figures 9(a) and 9(b). As it is seen, the NIM defect layer causes only one 
defect mode whose property is independent of the geometry of the structure. Moreover, the results 
show that the defect mode in the symmetric structure appears at a lower frequency compared 
to the asymmetric one. Furthermore, comparing the peak height of the defect modes in PIM 
and NIM defect layers shows that the peak height decreases when the NIM defect layer is considered. 



IV. CONCLUSION 

In this paper, the properties of defect modes coming from PIM and NIM defect layer at the center 
of a one-dimensional lossy MPC with symmetric and asymmetric structures on the transmission 
spectra have been investigated. Our numerical results show that defect modes with different 
characteristics for different geometries appear in the band gap. When the PIM defect layer is 
considered, there exists only one defect mode in the symmetric structure, but in the asymmetric 
structure two defect modes with different behaviors in the transmission spectra. Positions of the 
defect modes depend on the polarization of the waves and the incidence angles. As the incidence 
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angle increases the frequency of the defect modes moves toward higher frequencies regardless of 
the geometrical structure. Furthermore, one of the defect modes for TM polarized waves in the 
asymmetric structure disappears for incidence angles greater than . On the other hand, the results 
of the NIM defect layer show that only one defect mode which is independent of the geometric 
structure. Finally, an important result of this study is that by using a negative index material 
(metamaterial) in the structure of PCs, the defect modes show different behavior from the PCs 
with only positive index materials. In the latter conventional PCs case, as was mentioned earlier 
in the text, for a PIM defect layer, the symmetric geometric structure shows two defect modes 
and the asymmetric one shows one defect mode But in our case this is vice versa. Detailed 
analysis of the defect modes in metamaterial photonic crystals with symmetric and asymmetric 
structures will certainly provide useful information for designing new types of transmission narrow 
filters. 
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FIG. 1: Schematic diagrams of ID MPCs which arc embedded in air, with a defect layer, (a) asymmetric 
structure, (b) symmetric structure, where layer A is the NIM and layer B is PIM. The number of the periods 
is N and the incidence angle is 60 ■ 
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FIG. 2: Real parts of (a) the permittivity and (b) the permeability of layer A, e A and fx A , versus frequency 
for 7 = 0, 0.7, 1, and 1.5 GHz. 
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FIG. 3: Transmission spectra of TE polarized wave for the asymmetric ID MPC structures with PIM defect 
layer (solid line) and without defect layer (dotted line) for different incidence angles with 7 = 0.2 x 10~ 3 
GHz. 




FIG. 4: Transmission spectra of TM polarized wave for the asymmetric ID MPC structures with PIM defect 
layer (solid line) and without defect layer (dotted line) for different incidence angles with 7 = 0.2 x 10~ 3 
GHz. 
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FIG. 5: Frequency of the defect modes in the asymmetric ID PMC as a function of angle of incidence for 
both polarizations, with 7 = 0.2 x 10~ 3 GHz. 
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FIG. 6: Transmission spectra of TE polarized wave for the symmetric ID MPC structures with PIM defect 
layer (solid line) and without defect layer (dotted line) for different incidence angles with 7 — 0.2 x 10 -3 
GHz. 
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FIG. 7: Transmission spectra of TM polarized wave for the symmetric ID MPC structures with PIM defect 
layer (solid line) and without defect layer (dotted line) for different incidence angles with 7 = 0.2 x 10~ 3 
GHz. 
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FIG. 8: Frequency of the defect modes in the symmetric ID PMC as a function of angle of incidence for 
both polarizations, with 7 = 0.2 x 10~ 3 GHz. 
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FIG. 9: Transmission spectra of (a) asymmetric and (b) symmetric ID MPC structures with NIM defect 
layer (solid line) and without defect layer (dotted line) for normal incidence angle with 7 = 0.2 x 10~ 3 GHz. 
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